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Abstract. We present transverse momentum spectra, rapidity distribution and multiplicity of /1-hyperons 
measured with the HADES spectrometer in the reaction Ar(1.76A GeV)+KCl. The yield of E~ is calculated 
from our previously reported S~ /(A + £°) ratio and compared to other strange particle multiplicities. 
Employing a strangeness balance equation the multiplicities of the yet unmeasured £ ± hyperons can be 
estimated. Finally a statistical hadronization model is used to fit the yields of tv~ , K + , K°, K~ , <f>, A and 
E~ . The resulting chemical freeze-out temperature of T = (76 ± 2) MeV is compared to the measured 
slope parameters obtained from fits to the transverse mass distributions of the different particles. 
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1 Introduction 



Strange hadrons are particularly suitable probes of the 
high density phase of nuclear matter produced in few GeV 
heavy ion collisions. For instance, from systematic inves- 
tigations of subthreshold K + production tight constraints 
could be put on the nuclear equation of state at mat- 
ter densities of 2-3 po [JH1GIJH]. Furthermore, kaon phase 
space distributions and flow patterns are considered to 
be sensitive to the in- medium kaon potential [5H7]. On 
the other hand, due to strangeness conservation in the 
strong interaction, kaon production is intimately linked 
to the concurrent production of hyperons. While strange 
particle production is well understood in elementary NN 
collisions, in heavy ion reactions multi-step processes in- 
volving mesons or baryon resonances open up many ad- 
ditional production channels, even below threshold. Thus, 
strangeness-exchange channels like irA — > NK~ have been 
proposed to explain the observed K~ yields pi9"l[TfJ]. just 
as feeding through the — > K + K~ decay has been [TTj . 

Various aspects of strangeness production at SIS (Schw- 
erionen Synchrotron at GSI Darmstadt) energies have been 
investigated by the FOPI and KaoS experiments (for re- 
views see [T21HB] ). Evidently, any in-depth understanding 
of strangeness production and propagation in heavy ion re- 
actions requires information on all particles with open or 
hidden strangeness. The HADES collaboration has done a 
complete measurement in the system Ar+KCl at a bom- 
barding energy of 1.76 A GeV. Results on K° production 
have been published already in [13], on K + ,K~ and (f>- 
meson production in |11) and on the first observation at 
such a low beam energy of the double-strange E~ hy- 
peron in [TS]. To complete the picture, hyperon produc- 
tion remains to be addressed and this is the purpose of the 
present paper. We report here on the results obtained with 
the HADES detector on A production, from which, by ap- 
plication of strangeness conservation, we could estimate 
also the yield of the (not directly observed) S hyperons. 
Furthermore, we compare our set of particle yields to the 
result of a statistical hadronization model and discuss the 
implications. Note that the FOPI collaboration performed 
a similar analysis of strangeness production in the system 
Ni+Ni at 1.93A GeV [16]. 

In Section 2 of this paper we give first a brief overview 
of the HADES detector and relevant details of the Ar+KCl 
data taking and then proceed to describe the employed 
particle identification and A reconstruction procedures. In 
section 3 we present spectra and production yields of the 
A hyperons. In section 4 the A result is used to extract 
the yield of the double-strange £~ from our previously 
published S~ / A ratio [TS]. With all experimental yields 
established, strangeness balance is applied to estimate the 
yield of the unobserved charged S hyperons. We compare 
all yields obtained in Ar+KCl with respect to a statisti- 
cal hadronization model and confront the fitted chemical 
freeze-out temperature with the measured slope parame- 
ters obtained from transverse mass distributions. Finally 
we summarize our findings in section 5. 




impact parameter [fm] 

Fig. 1. Impact parameter distributions of all and LVL1 se- 
lected Ar+KCl reactions obtained from the UrQMD transport 
code [IE]. 

2 Experimental setup 

HADES is a charged-particle detector consisting of a 6- 
coil toroidal magnet centered on the beam axis and six 
identical detection sections located between the coils and 
covering polar angles between 18° and 85°. Each sector is 
equipped with a Ring-Imaging Cherenkov (RICH) detec- 
tor followed by Multi-wire Drift Chambers (MDCs), two 
in front of and two behind the magnetic field, as well as a 
scintillator hodoscope (TOF/TOFino). Lepton identifica- 
tion is provided mostly by the RICH and supplemented at 
low polar angles with Pre-SHOWER chambers, mounted 
at the back of the apparatus. Hadron identification, how- 
ever, is based on the time-of-flight and on the energy-loss 
information from TOF/TOFino, as well as from the MDC 
tracking chambers. A detailed description of HADES is 
given in [17] . 

An argon beam of ~ 10 6 particles/s was incident with 
a beam energy of 1.76 A GeV on a four- fold segmented 
KC1 target with a total thickness corresponding to 3.3 % 
interaction probability. A fast diamond start detector lo- 
cated upstream of the target was intercepting the beam 
and was used to determine the time-zero information. The 
data readout was started by a first-level trigger (LVL1) re- 
quiring a charged-particle multiplicity, MUL > 16, in the 
TOF/TOFino detectors. Based on a full GEANT simula- 
tion of the detector response to Ar+KCl events generated 
with the UrQMD transport model [T5] , we found that the 
event ensemble selected by this (LVL1) trigger condition 
has a mean number of participating nucleons {(A part }) 
equal to 38.5 + 3.9. Figure [T] illustrates the impact param- 
eter distributions obtained from UrQMD calculations for 
two event selections: all inelastic events and according to 
the experimental LVL1 trigger condition. 

The particle identification was done by a velocity vs. 
momentum x polarity correlation, where the velocity was 
determined by the time-of-flight measurement in the TOF 
and TOFino scintillators with respect to the time-zero in- 
formation and the tracked flight path. If needed, addi- 
tional particle discrimination was gained from the energy- 
loss (dE/dx) information in the MDC and scintillators. 
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3 A hyperon yield and spectra 
3.1 A identification 

The particle identification of kaons, 0, tt~ and E~ is de- 
scribed in [TTlfTMn] . Here we add only those details spe- 
cific to the reconstruction of A hyperons in their decay 
channel A -> p + n~ (B.R. = 63.9%, cr=7.89 cm [19]). 
Note that at our beam energies the reconstructed A yield 
contains also a contribution from decays of the slightly 
heavier £° baryon into a A and a photon. 

The decay products of the A hyperons have been iden- 
tified using the MDC dE/dx and time-of-flight informa- 
tion. The topology of the A decay into p-7r _ pairs has 
been used to suppress the combinatorial background of 
uncorrelated pairs. Cuts on the distance between the pri- 
mary event vertex and the decay vertex (dyo), on the dis- 
tances between the proton (d p ), respectively the it" (d^- ) 
track and the primary vertex, on the distance of closest 
approach between the two tracks (ddca) and on the dis- 
tance of the reconstructed mother particle trajectory to 
the primary vertex (d p7r - ) were applied. Furthermore, a 
minimum opening angle (a pir -) was required to guaran- 
tee a good decay vertex resolution. All selections used for 
the analysis are listed in Table [TJ For an estimation of the 
systematic errors, the cut values have been varied within 
reasonable limits. In total, 28 different cut combinations 
were thus investigated, resulting in a total amount of re- 
constructed A hyperons ranging from 36k to 191k. The 
full reconstruction chain with all corresponding efficiency 
corrections was applied for each of the cut variations. All 
systematic errors are based on these cut variations and 
represent the maximum/minimum deviation to the results 
obtained with the chosen cut values in Table Q~] 

Figure[5]shows the invariant-mass spectrum of all proton- 
7r~ pairs which passed the cuts listed in TableQ] An event- 
mixing technique has been used to model the combinato- 
rial background of uncorrelated pairs. Displayed in Fig. [2] 
as a grey shaded histogram, the mixed event background 
was normalized to the data on the left and right side of the 
A peak. In total, for the optimal cut selection listed in Ta- 
ble[TJ about 100000 A hyperons were reconstructed, with a 
mean signal-to-background ratio of 0.56. From a Gaussian 
fit to the peak, the pole mass is determined to be 1114.3 
MeV/c 2 , i.e. about 1.4 MeV/c 2 away from its listed value 
[Hi] . We attribute this small difference to residual deficien- 
cies of our track reconstruction and detector alignment. 



3.2 A spectra 

For further kinematical studies the A signal has been de- 
termined in nine rapidity bins, ranging from —0.75 < 
Vc.m. < +0.15 in steps of 0.1, and up to ten transverse 
mass bins in steps of 50 MeV/c 2 . The background sub- 
tracted signal yields were corrected for acceptance and re- 
construction efficiency using a full GEANT simulation of 
the detector system described in [14] and a track-embedding 
method. The geometrical acceptance, which also includes 
the branching ratio of A — > p + ir~ of 0.639, shows a 
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Fig. 2. Top: Invariant mass of all identified proton and n~ 
pairs after several cuts on the topology of the A decay kine- 
matics were applied (see text for details). The grey shaded his- 
togram shows the mixed-event combinatorial background, nor- 
malized to the signal spectrum between 1080-1100 and 1130- 
1150 MeV/c 2 . Bottom: A signal after background subtraction; 
the solid red line shows a Gaussian fit to the signal. 



smooth behaviour as a function of the transverse mass 
and varies for most of the bins between 13% and 34%. It 
is defined by the requirement that both daughter parti- 
cles have hits in all MDC planes. The A reconstruction 
efficiency is composed of the single track reconstruction 
and particle identification efficiencies (~80% per track) 
and the cuts on the yl-hyperon decay topology. The latter 
one clearly dominates the reconstruction efficiency which 
has values of 3% to 10%. The dominant topology cut is the 
one applied on the distance between the primary vertex 
and the /1-decay vertex. It is in the order of the /t-hyperon 
mean decay length (see Table QJ. Acceptance and recon- 
struction efficiency are plotted for the mid-rapidity region 
in Fig. H 

The acceptance- and efficiency-corrected transverse- 
mass spectra of A for the various rapidity bins are pre- 
sented in Fig. 21 Shown is the number of counts per LVL1 
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Table 1. Topological conditions values chosen for the A analysis (see text). 
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Fig. 3. (Color online) Acceptance and reconstruction efficiency 
at mid-rapidity for yl-hyperons as a function of the reduced 
mass. The plotted acceptance (red open circles) includes al- 
ready the branching ratio of 0.639 of A — > p + n~ . The recon- 
struction efficiency (black filled circles) consists of the single 
track reconstruction efficiencies of the daughter particles and 
the efficiencies of the topological cuts applied to improve the 
signal-to-background ratio (see text for details). 



trigger, per transverse mass and per rapidity bin, divided 
by m 2 . This representation is chosen in order to easily 
apply Boltzmann fits to the resulting distributions, ac- 
cording to 
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The solid lines in Fig. [4] show the results of Boltz- 
mann fits, where TB{y c .m.) represents the inverse slope 
of each distribution. The resulting TB(y c .m.) values are 
then plotted in Fig. [5] as a function of the center-of-mass 
(cm) rapidity (y c . m . = V — y{cm), where y{cm) = 0.858 
for symmetric collisions at 1.76 A GeV). The full symbols 
display the measured data, whereas the open ones are the 
data reflected at cm. rapidity. The error bars represent 
the statistical errors. Assuming a thermal source, these 
temperatures are expected to follow the relation 
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yielding an effective temperature of T e ff — (95.5 ±0.7(stat.) 
+2.2(syst.)) MeV. Here the systematic error corresponds 
to the variation of the cut values described above. 

For each rapidity bin the transverse mass spectrum 
was integrated in the following way: the yields in the cov- 
ered bins were added and the fits were used to extrapolate 



10 
10 



-2 



-3 



1Z io-' 

CM 

O 

=J 10 
a> 

^ 10 



-6 



-7 



TJ 



10 

10"' 



10 



^ 10" 
■o 

^ 10 

CM 4~ 

E 



10 



■11 



10" 12 
10" 13 
10" 14 
10" 15 




A +0.05<y cm <+0.15 (x10 )Q -0.45<y c m <-0.35 (x10 ) 
ft -0.05<y°J<+0.05 (x10 7 ) O "0-55<y^ <-0.45 (xloT 
: A -0.15<y°%-0.05 (x10 6 ) )K -0.65<y^ <-0.55 (x10 1 ) 
Fa -0.25<y' %"0-15 (x10 5 ) V -0.75<y c m <-0.65 (x10T 
: -0.35<y^<-0-25 (x10 4 ) ^ 



200 

m t -m A 



400 
[MeV/c 2 ] 



Fig. 4. Reduced (m t — tua) transverse mass spectra for dif- 
ferent rapidity selections. For better legibility, the spectra are 
scaled as indicated in the legend. The solid lines are fits with 
Eq.[T]to the data. 



into the unmeasured kinematic regions. The fits were in- 
tegrated from to the lower bin edge of the first measured 
point and from the upper bin edge of the last measured 
point to infinity. The fraction of the extrapolated yield in 
the transverse mass spectra to the total yield is 36-43% for 
the rapidity bins in the range —0.65 < y c .m. < 0.15, and 
65% for the rapidity interval (—0.75 < y c . m . < —0.65). 
The results are shown in Fig. |6l where the obtained ra- 
pidity density distribution is displayed. The full trian- 
gles show the values calculated by the integration of the 
transverse mass spectra, while the open triangles repre- 
sent points reflected with respect to the center-of-mass 
rapidity. 

For the determination of the total A multiplicity per 
LVL1 event, the measured spectra were integrated. The 
extrapolation into the unmeasured region was done by fit- 
ting either a gaussian or a linear function to the first four 
data points, as shown in Fig. [S] The mean value of these 
two different extrapolations is used for the total yield. 
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A detailed description of the A analysis can be found in 

EDI. 



4 Discussion 
4.1 Particle yields 

Knowing the yield of the A+E° hyperons and the S~ /(A+ 
E°) ratio yields the production rate of the E~ to (2.3 ± 
0.9) x 10~ 4 , adding statistical and systematic errors quadrat- 
ically. Note that this value is of the same order of magni- 
tude as the yield of the (j> meson [TTj . 

Table [2] summarizes all particle yields extrapolated to 
full phase space as well as the corresponding inverse slope 
parameters from fits to the particle m t spectra; results on 
K+,K-,K%, and <t> are taken from [TTJCE1]. 



cm. 



Fig. 5. Inverse-slope parameters from fits with Eq. |T} to the A 
hyperon transverse mass spectra as a function of rapidity. The 
solid line is a fit with Eq. © to the data points. The effective 
temperature T e / / is the function value at mid rapidity. 



x10 



-3 



I 




A A (measured) 


A : 


: a A (reflected) 




A A A A 








L I \ 
I \ 




* 

i 

/ i 

f' ; 


i : - 


-1 


1 


^c.m. 





4.2 Strangeness balance 

The strong interaction conserves strangeness, i.e. the num- 
bers of s and s quarks produced in a heavy ion reaction 
must be equal. As those quarks are ultimately bound in 
hadrons the multiplicities of strange particles fulfill a bal- 
ance equation which can be written at SIS energies as: 



50 
40 



■o 30 



20 
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Fig. 6. Rapidity distribution of A hyperons. The closed sym- 
bols refer to measured data points calculated from the trans- 
verse mass spectra, whereas the open symbols show the data 
points reflected about center-of-mass rapidity. For the extrap- 
olation to unmeasured rapidity values a linear and a Gaussian 
function were fitted to the first four data points (see text). 



The fraction of the Gaussian extrapolation to the total 
yield is about 4.2%, whereas the fraction of the linear ex- 
trapolation is negligible. The inclusive total A multiplic- 
ity per LVL1 event was found to be (4.09 ± 0.1(stat.) ± 
0.17(extr.)^Q37(syst.))xl0~ 2 , where the second error refers 
to the extrapolation uncertainty in m t and the third one 
to the systematic error obtained from the cut variations. 



K 



£ 0± + A + K~ + K° + 2H ' 



(3) 



where, for simplicity, the symbols denoting the particles 
stand for their respective yields at the time of production. 
Note that this equation takes care of the strong decay of 
heavier strange resonances via the counting of their de- 
cay products, namely kaons and As. As mentioned ear- 
lier, the U° can not be separated from the A, thus this 
contribution is to be counted explicitly together with As. 
Analogously, according to our analysis procedure, most 
decay products feed the A channel and are counted 
as As, therefore the factor in Eq. © in front of the 
is two instead of four. (Note that anyhow the contri- 
butions are small.) In case of the neutral kaons, we mea- 
sure in fact the yield of the K®, which obeys the equality 
= (K° + K°)/2. Assuming isospin symmetry the yield 
of the K° should be contributing here at the same order 
as the K~ yield. Eq. can then be rewritten using the 
measured yields. Hence the unobserved hyperon yield 
can be estimated as: 



2K° S - {S° +A) — 2~ - - 3if~ (4) 



Still heavier multi-strange particles, e.g. Q hyperons, have 
significantly higher production thresholds and should not 
contribute sizeably at SIS energies. From Eq. (HJ a total 
multiplicity of charged £ hyperons of (7.5 ± 6.5) x 10 -3 is 
deduced when using the values of the multiplicities listed 
in Tabled The error is the quadratic sum of the statistical 
and systematic error of the different yields. If one assumes 
isospin symmetry for the U^'° yields one can subtract 
the E° contribution from the A yield and finds the ratio 
Z ,± '°/yl is 0.3 ± 0.26 although the difference in mass is 
only 10%. 
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Table 2. Multiplicities (i.e. yield/LVLl event) and effective temperatures of particles produced in Ar+KCl reactions at 1.76A 
GeV. The error on the E and S~ yield is the quadratically added statistical and systematic error. 



Particle 


Multiplicity 


Teff [MeV] 


Reference 




3.9 ±0.1 ±0.1 


82.4±0.1+^g 


m 


A + E" 


(4.09 ± 0.1 ± 0.17^0.37) x 10~ 2 


95.5 ±0.7 ±2.2 


this work 


K+ 


(2. 8 ±0.2 ±0.1 ±0.1) x 10~* 


89 ± 1 ± 2 


11 


K°s 


(1.15 ±0.05 ±0.09) x 10~ 2 


92 ± 2 


14 


K~ 


(7.1 ± 1.5 ±0.3 ±0.1) x 10~ 4 


69 ± 2 ± 4 


ED 


<P 


(2.6 ±0.7 ±0.1 -0.3) x 10~ 4 


84 ±8 


ED 




(2.3 ±0.9) x 10~ 4 




.15. 


E+ + E~ 


(0.75 ±0.65) x lCT* 




estimated via strangeness balance 



The only other published multiplicity of charged S hy- 
perons in heavy ion collisions, based on a similar analysis 
of strangeness yields measured with the FOPI detector at 
GSI in Ni±Ni reactions at 1.93A GeV, is (7± 8 + 32 -17) x 
10~ 3 [16]. Differences with respect to our analysis are: (1) 
a higher beam energy (1.93 vs. 1.76A GeV), (2) a larger 
reaction system (58+58 vs. 40+37), and (3) a different 
centrality selection ((A part ) — 71 vs. 38.5) resulting in a 
higher A yield (0.137 ± 0.005 + 0.007 - 0.008) and inverse 
slope parameter (119±l+9— 7). In view of the higher bom- 
barding energy and larger system size, one would expect a 
larger charged S contribution. The resulting Z ,±, ° / 'A ra- 
tio is 0.08 ± 0.09 + 0.33 - 0.18. Unfortunately, the large 
uncertainties prevent us to draw a firm conclusion on the 
behavior of the respective cross sections with energy. 

4.3 Comparison with statistical hadronization 

Statistical hadronization models (SHM) have been suc- 
cessful in fitting particle yields or yield ratios from rel- 
ativistic and ultrarelativistic heavy ion collisions [21,22, 
UllMj- With the help of SHM fits it has been possible 
to reconstruct systematically the chemical freeze-out line 
in the T — fib plane of the nuclear phase diagram with 
Hb being the baryochemical potential (see e.g. [2"4"ll2"o] ) . 
However, while the various SHM approaches agree fairly 
well at high bombarding energies, discrepancies appear 
in the low-energy regime. Indeed, at the lower energies 
it is not even clear, whether chemical equilibrium can be 
reached [57] and therefore the question arises whether a 
statistical treatment of particle production is meaning- 
ful. The situation is further complicated by the need for 
strangeness suppression, which is handled differently in 
the various SHM implementations. Furthermore, at SIS 
energies, only pions are produced abundantly. Heavier and 
especially strange particles are rare, and their yields were 
mostly poorly known. Hence in the past only few particle 
yields with small statistical errors were available as input 
to the fit procedure. In the following we fit eight particle 
yields obtained from our Ar+KCl run with a statistical 
hadronization model. 

We choose the freely available THERMUS code [28], 
using the mixed canonical ensemble where strangeness is 
exactly conserved while all other quantum numbers are 



calculated grand canonically. We handle the strangeness 
suppression by introducing a strangeness correlation ra- 
dius R c within which strangeness has to be exactly con- 
served; this is discussed in (29j . We fit simultaneously all 
particle yields listed in Table [2] except for the U^, as 
well as the mean number of participants (A part ) and con- 
strain the charge chemical potential fiQ using the ratio 
of the baryon and charge numbers of the collision sys- 
tem. We find the chemical freeze-out at a temperature of 
Tchem = (76 ± 2) MeV and at a baryochemical poten- 
tial of Hb = (799 ± 22) MeV. The strangeness correlation 
radius comes out as R c = (2.2 ± 0.2) fm, which corre- 
sponds to about half the fitted radius R = (4.1 ± 0.5) 
fm of the whole fireball. The exclusion of the E~ from 
the fit changes the parameters only on the percent level, 
but the x 2 /d.o.f. value of the fit improves from 13.9/4 to 
7.8/3. Fig. 0shows the resulting freeze-out point together 
with a compilation of similar points [30,24,31 in the T - 
Hb plane. Our result, as well as the FOPI result from the 
collision system Al+Al at 1.9A GeV, differ from the reg- 
ularity of freeze-out points following the fixed energy per 
particle condition (E) / (N) w 1 GeV, which is one of the 
commonly proposed freeze-out criteria [24] . This might be 
due to the light collision systems, since small systems have 
the tendency to show higher freeze-out temperatures |25) . 

A detailed comparison of the data with the statistical 
model fit is shown in the upper part of Fig. [S] while the 
lower part of this figure depicts the ratio of data and fit. 
All particles except for the E~ are well described. 

A particularly interesting case is the <f> meson. The <fi 
is treated as a strangeness neutral object in the R c for- 
malism and is therefore not suppressed at all. Its yield is 
well described by the SHM. This means that the 4> yield is 
compatible with the assumption that it takes part in the 
equilibration of the hadrons. This is quite different from 
the situation at higher bombarding energies, where the cj> 
requires indeed an effective strangeness between 1 and 2 
to have the appropriate suppression in the SHM and to 
reproduce the data [29] . For an understanding of <p pro- 
duction, one may have a look at the 4>/K~ ratio which, 
according to the SHM with R c , should rise at low beam 
energy. Such a behavior is indeed supported by our data, 
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Fig. 7. (Color online) Chemical freeze-out points in the T - 
fib plane. The filled black circles (a) are taken from [30], the 
black open triangles (b) are from [24]. The red circle is taken 
from [3T]. The THERMUS fit to our Ar+KCl data is shown 
as blue triangle. The dashed line correponds to a fixed energy 
per nucleon of 1 GeV, calculated according to [24] . 



10 2 
10 



210" 1 
o 

10" 3 



,-5 



10 



§ 1.5 

m 

w 1 



• Data, Vs NN =2.61 GeV 

T=75.8±2.1 MeV, |i =799.4+21 .6MeV 
R„=2.2±0.2fm, R =4.1±0.5fm 



24 + 9 



A part 7i- A K + K $ ST E + " 



Fig. 8. (Color online) The upper plot shows the yields of sec- 
ondary hadrons in Ar+KCl reactions (filled red circles) and 
the corresponding THERMUS fit (blue bars). The lower plot 
shows the ratio of the experimental value and the SHM value. 
For the E~ the ratio number is quoted instead of a point. 



as already discussed in (11) . and the ratio seems to ap- 
proach the value seen in elementary NN reactions 32. . 

According to the strangeness suppression mechanism 
implemented in SHM, the double-strange S~ (S=2) should 
be suppressed strongly with respect to the <p with its hid- 
den strangeness (S=0). Nevertheless, our measured S~ 
yield is of the same magnitude as the one of the 0, i.e. the 
data show no indication for any strangeness suppression. 
This is very surprising since the S~ yields observed above 
threshold at RHIC [33], at SPS [3J] and even at AGS [35] 
are consistent with statistical model fits. In fact the same 
secondary pion-hyperon process 7r + Y — > (f> + Y, which 
was invoked by Kolomeitsev and Tomasik |36] to explain 
the enhanced 4> yield, can here be the origin of the high 3 
production via the reaction tt + Y — > 3 + K . To get a bet- 
ter understanding, we may have to move away from the 
SHM. One may calculate the probability for the produc- 
tion of two ss pairs in one collision. Assuming that both 
pairs are independently created, their production proba- 
bility P2ss is given as the square of the single-pair pro- 
duction probability P s -g. Keeping associated production 
in mind, P S s can be estimated as the combined multi- 
plicity of all particles that carry a strange quark, respec- 
tively the combined multiplicity of all anti-strange par- 
ticles, i.e. K + + K° + (f), yielding P S g ~ 0.05 and hence 
P 2s s — 0.0025. Considering that the observed 3~ yield 
is in fact an order of magnitude smaller, we conclude 
that in 10% of these events both s quarks end up to- 
gether in a 3~ , whereas from strangeness suppression in 
the SHM one obtains less than 1%. A different realiza- 
tion of the SHM using the strangeness canonical ensem- 
ble and 7 S for additional strangeness suppression delivers 
comparable freeze-out parameters, with T c i iem = (76 ± 5) 
MeV, /i b = (791 ± 33) MeV, R = (4.1 ± 0.9) fm and 
7 S = 0.37 ± 0.04 but fails to reproduce the <fi multiplicity 
by an order of magnitude due to its suppression with 7^. 



4.4 Chemical vs. kinetic freeze-out 

The temperature T c h em obtained for the chemical freeze- 
out can be compared with the inverse-slope parameter 
T e f f obtained from Boltzmann fits to the rrit spectra of the 
different particle species. Apparently most of the inverse- 
slope parameters are higher than the chemical freeze-out 
temperature of the system. A pure Boltzmann shape can 
be distorted by various effects, like collective motion or 
early vs. late particle decays. One example is apparent in 
the difference between the K + /K® and K~ slopes (see 
Fig. [9]). The much lower value of T e ff of the K~ has often 
been interpreted as due to its much later freeze-out time 
j!2l neglecting the admixture of soft K~ stemming from 
(j> decays. In [37] however it was shown that these soft K~ 
indeed affect the shape of the spectra. 

Effects of collective flow, on the other hand, should 
influence the transverse mass slope more, the higher the 
particle mass. From Fig. [9] where the fitted temperatures 
are ordered by increasing particle mass, this seems not to 
be a strong effect as expected for a small collision system 
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Fig. 9. (Color online) Effective temperature T e // of all mea- 
sured particle species as a function of their mass. The horizon- 
tal line and error band show the chemical freeze-out tempera- 
ture Tchem from the THERMUS fit, whereas the dashed-dotted 
line and the red error band show a linear fit to the data points 
(K~ are excluded, see text for details). 



like Ar+KCl. However the inverse-slope parameter seems 
to be slightly decreasing with decreasing mass. 

To take this effect into account a linear fit to the data 
points was applied (dashed-dotted line and red error band 
in Fig. in]). The K~ were excluded from the fit for the rea- 
sons discussed above. The fit is clearly dominated by the 
A and K + /K® data points with small relative errors com- 
pared to the <p and tt~ data points. Within errors, the 
extrapolated fit value at mass is still above the chemical 
freeze-out temperature. Hence the presented data implies 
an inversion of the kinetic- and chemical freeze-out sce- 
nario, which cannot be the case for obvious reasons. This 
means either the statistical model approach for small reac- 
tion systems and small energies is not applicable, and/or 
the unique kinetic freeze-out for all particles with one ra- 
dial flow velocity is a too naive assumption for this system. 
Future measurements of HADES in the reaction systems 
Au+Au and Ag+Ag will give more insight into the com- 
plex dynamics at low energies. 



5 Summary and conclusions 

We have presented phase space distributions of A hyper- 
ons in Ar+KCl at 1.76 A GeV measured with the HADES 
spectrometer at GSI. Combining the measured A + 27° 
yield with our former data on strangeness production in 
this system we have estimated the yield of the double- 
strange E~ hyperon. We find that it is of the same order of 
magnitude as the one of the (j> meson. The fraction of the 
unobserved charged 27^ hyperons could be constructed 
using strangeness conservation. 

Applying a statistical model fit to these hadron yields, 
a fair agreement, except for the S* - , in a strangeness- 
canonical approach is achieved. The </>, however, is well 
reproduced without any suppression, in sharp contrast to 



the situation at higher energies, where a suppression is 
observed. 
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